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The effect of light quark masses, as seen from short distances, is discussed. The large logs of m q /A are shown to simulate 
"vacuum condensates". 



1. The aim of this Letter is to show that the role of light 
quark masses has heretofore been strongly underrated. Indeed, 
according to the popular philosophy the effect of light quark 
masses should be relatively unimportant, while the leading role 
is ascribed to the complex non-perturbative vacuum which is 
considered as only lightly affected by m q J=0 [1]. Nobody, 
however, has as yet seriously studied the effect of large logs of 
m q l fi, fl being the renormalization parameter, which appear 
in perturbative calculations. Below we argue that these logs, 
when resummed, can drastically change the naive understand- 
ing of the limit m q —> , so that the latter should rather be 
considered in the spirit of Bogoliubov's quasi-averages [2]. By 
this we mean that the real QCD near m q ~ may differ quali- 
tatively from the theory with massless quarks if such a theory 
exists (which is far from being obvious). 

To study the limit m q — » of the full theory one must first 
expand the full Green functions in m q and only then try, if one 
wishes, to calculate the resulting coefficients by expanding 
them in a s ■ For obvious reasons such a procedure is hardly 
possible. Instead one usually proceeds in a completely opposite 
direction: one takes the limit m q — > , graph-wise, after ex- 
panding in a s , so that \og(m q l ft) 's are always suppressed by 
powers of m q . 

It would be preferable to try to see the whole truth by re- 
summing, as best we could, those logs prior to making conclu- 
sions. Needless so say, this is no substitute for a non- 
perturbative analysis, but it is reasonable to expect that this 
way corresponds better to the original intention. 

2. The Bogoliubov-inspired [3] extension principle (EP) for 
studying asymptotics of Feynman integrals [4] in conjunction 
with the R* -operation for removing a class of infrared diver- 
gences in the presence of ultraviolet ones [5], offer a means to 
study the structure of the m q — > expansion of Feynman inte- 
grals straightforwardly and to all orders of perturbation theory 
(PT). 

To avoid irrelevant combinatorial complications we take a 
very simple example, but which retains all the essential fea- 
tures of the full proof 1 . Our example is similar to the one con- 
sidered in [4] and our reasoning closely follows [6]. Therefore 
only formal steps of the derivation will be indicated below. 
The MS-scheme [7] is used throughout the paper. 

Consider a free theory with a massless field (p and a mas- 
sive one $ in 4 dimensions, and the current 
j(x) = :(p(x)<S>(x): . Consider 



U(Q\m 2 )=ije^d D x(RTj(x)j(0)) Q 
d D p 1 



= 1 



{2k) d (p-Q) 2 (m 2 +p 2 ) 



■C.T. 



(1) 



where Q 2 = -q 2 , and the rotation to Euclidean region has 
been performed on the r.h.s. of (1). The UV counterterm in (1) 
does not depend on m . Formal expansion gives: 



Eq.(l) = J 



d D p 1 
{2n) D (p-Q) 2 



C.T. + o(m 2 ) . 



(2) 



The UV counterterm still subtracts the UV divergences from 
the first term in (2). The second term has no UV divergences 
but a logarithmic IR one. Following EP [4], to obtain a correct 2 
expansion in m, we must add a counterterm into (2): 



~^~+z{m){2n) D 8{p). 
P P 



(3) 



(NB Derivatives of 5(p) need not and should not be added 
[4], [6].) To fix the unknown z(m) we demand that 



(2w) 



1 



2 2 

m +p 



1 m l 
P P 



■z(#n) (2x) D 8(p) 



(4) 



for any smooth \ff(p) such that \jf(Q) = 1 and lim \jf{p) < +°° 



*0 



(5) 



[4], [6]. Choosing \jf(p) = 1 we obtain: 

z{m) = \ d ° p n * 2 ^(r<£ 2 (0)) 

J (2n f m 2 +p 2 \ lo 
(note the absence of normal ordering in <I> 2 (0) ). 

Using (RT<3> 2 (0)) = n D - 4 m 2 Z + (T<£> 2 (0)) to express 
(5) via convergent quantity (RT<i> 2 (0)) we finally obtain: 

n(G,m) = n(G,o) 

m 2 C 1 (Q,n) + (Rm 2 (0)) C 2 (Q,n) 



Q 



- + o(m z ). 



(6) 



Each quantity in (6) is finite by construction. 
Interaction modifies (6) in three respects: 

(i) C\ and C 2 become series in coupling constants and 
\og(Qlft) ; 

(ii) (RT<t> 2 (0)) Q goes over into {^ 2 (0)}=(RT[<P 2 (0)S]) Q *0 , 
S being the S-matrix; 

(iii) along with {$2 (Q)} , there appear the vacuum averages 
{(p 2 (0)} and {(p(0)<I>(0)} with their own coefficient functions, 
due to a mixing mechanism which is analogous to the mixing 
of operators under UV renormalization. 

A generalization of (6) reads: 

i\e i " x Ax{A{x)B{Q)} = £ £ C^ft) m 2 '{o (0)}. (7) 

e 2 ^+~ i=0 j L 1 

{Oj(0)} are vacuum averages of local composite MS- 
renormalized operators [8], inevitably non-vanishing at least 
while m * . The leading power behavior of Qj at Q — » °o is 
determined by power counting and is modified by jj. only via 
the log(g/M) contributions. For explicit formulae for calcu- 
lating Cjj(Q,n) see [9]. 



1 The full and general proof will be given elsewhere; it repeats, in essence, 
the derivation below step by step. 



2 "Correct" means (i) "IR convergent" and (ii) "ensuring the o(m 2 ) esti- 
mate for the neglected terms in (2)". 
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3. Several conclusions can be drawn from the above deriva- 
tion: 

(i) The techniques of dimensional regularization and mini- 
mal subtraction [7] are crucial for achieving the neat separation 
of Q and jl in (7) (the phenomenon of "perfect factorization" 
in terms of [6]; see also [10]). 

(ii) The quantities {Oj} in (7) summarize information on 
the small-p behavior of the theory. So the above derivation is, 
in fact, a formalization of heuristic arguments of [1] which led 
to introduction of non-zero vacuum condensates into deep- 
Euclidean QCD amplitudes (for a renormalization group analy- 
sis of {Oj} see below). 

(iii) By re-expressing the MS-renormalized operators in (7) 
in terms of, say, Zimmermann's normal products [1 1] with zero 
vacuum averages, all the \og(Q//l) terms could only be 
pumped into the coefficient functions but never got rid of. It 
should be clearly understood that the vacuum condensates 
{Oj} are factorization-prescription dependent quantities, and 
the virtue of the MS-scheme consists not in making them non- 
zero, but in collecting within {Oj} all the logs of mill. 

(iv) It is clear from the derivation that if there are heavy par- 
ticles in the theory and/or there are other large momenta (in the 
case of three- etc. -point functions), all the heavy, i.e. 0(q 2 ) , 
parameters will enter only coefficient functions in (7), so that 
(7) represents a general form of how the world of "light" con- 
fined particles affects the "heavy" world. That is to say, when- 
ever the introduction of vacuum condensates is justified from 
the phenomenological viewpoint of [1], they do appear as an 
effect of light quark masses already at the perturbative level. 

4. Now, the question arises: does the perturbative analysis 
of the vacuum condensates {Oj} make sense? 

Consider the set O d of the aggregates m 2i {Oj} with a 
given total dimensionality d in units of mass (they contribute 
to the same power of Qr 1 and form a closed set with respect 
to renormalization). The RG equation for O d reads: 

d-y{a s ) 



3 x 3 



o d = o, 



where y is the matrix of anomalous 
dm 



d jU 



(8) 

dimensions, 
(9) 



m B ,a SB 

Eq.(8) is similar in form to the well-known evolution equa- 
tions for moments of structure functions of DIS (for a review 
see e.g. [12]). The solution to (8) reads: 



O d (m) = O d (mA — 
where a*(m) is 



a (m) 



a* (mo) 



x[l + a 2 a*(m) + ...], (10) 
/}* • Assuming that 



governed by 

a*(m) « (fclog(m 2 /A* 2 )) _1 with A*~Ams -100-^500 MeV , the 
conclusion is unavoidable that with O d (m) and light quarks 
we are in a deep non-perturbative region. 

5. The most important point perhaps is that even if the m d 
factor in (10) prevails in the chiral (i.e. m q — > 0) limit, e.g. if 
(qq) in the above interpretation vanishes as m u d — > , such 

a possibility is by no means inconsistent with the current- 
algebra-PCAC phenomenology, as is explained in the appendix 
E of [13]. As explained there, if (qq) Q vanishes in the chiral 



limit then the Gell-Mann-Okubo formula for the r\ -meson is 
lost, while the masses of the quarks u and d are larger than in 
the standard framework, which gives one a possibility to ac- 
count for the large deviations from the chiral-invariant picture 
observed in the pion-nucleon scattering. 

However, theoretical estimates of the vacuum condensate 
values from (10) are hardly feasible at present — one of the 
reasons is that the value of A /A* is related to the strong- 
coupling behavior of fi m , which is unknown. 

6. To conclude, we have demonstrated that upon resumming 
large logs of the ratio m q /fl , where m q represents light quark 
masses and jl is the renormalization [= factorization] parame- 
ter, within usual perturbation theory, a collection of "vacuum 
condensates", i.e. non-zero vacuum averages of composite op- 
erators, emerges whenever they are liable to appear according 
to, and in exactly the same form as prescribed by, the phe- 
nomenological rules of [1]. The renormalization group analysis 
shows that with those condensates we are in a strong-coupling 
region, so that their behavior at m q — » is hardly predictable 
at present. However, the fact that the \og(m q / fl) 's fit into the 
pattern of vacuum condensates is in itself very nice, for it en- 
sures that those uncalculable terms, being hidden within con- 
densates, cannot modify the formulae of [1], whatever is the 
interpretation. 
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Comments (December, 1998) 

The subject of the comments below is to explain the sys- 
tematic approach to studying the structure of power corrections 
which was demonstrated in this Letter and which remains the 
only valid and completely meaningful scenario for studying the 
structure (not necessarily numerical values) of power- 
suppressed corrections. 4 

Bibliographic comments 

This Letter was in the first burst of publications ([4], [6], 
[9]) on what eventually evolved into the theory of Asymptotic 
Operation (hereafter AO; see [19] and refs. therein, and [25] 
for the more recent ultimate extension to arbitrary non- 
Euclidean situations). The formulas discovered within the 
framework of this vast and powerful theory formed the basis of 
an entire industry of multiloop calculations (NNL moments of 
DIS structure functions, [23] and many subsequent publica- 
tions; mass expansions, [30] and refs. therein 5 ). 

Remarks on the general method 

The derivation presented in sec. 2 follows the pattern of AO 
explicitly worked out to all orders in [14], [15] (for systematic 
discussions see [15] and the review [19]). At the time of writ- 
ing of this Letter, the notion of AO had not yet crystallized out 
from concrete calculations, instead I emphasized the "extension 
principle" which is a central but somewhat abstract element in 
prescriptions for writing out formulas of AO for concrete cases. 
As to the R* -operation mentioned in the first paragraph of 
sec. 2, it was actually obtained with an essential use of the 
method of AO (as explained in [19]) but final formulas could 



4 Concerning the discussion of the chiral limit in sec. 5 of the main text 
I only add that, given how nearly impossible it is to do and publish work 
that goes against a prevailing dogma created by many fine experts, I would 
not be surprised if the possibility of quark condensates vanishing in the 
chiral limit (although at a slower than canonical rate; cf. (10)) would turn 
out to be correct in the end. But that is only a possibility which one can 
neither prove nor completely exclude. 

5 Unfortunately, the publications by Kiihn et al. systematically fail to give 
proper credit to the original publications in which the mass expansion for- 
mulas they used were discovered, namely, [14]— [1 8]. This is the more de- 
plorable that the resulting misinformation of the research community as 
regards the true source of the underlying mathematical expertise proves to 
be detrimental to the progress of theoretical high energy physics because 
the full potential of the theory of AO is far from being exhausted (cf. [31]). 



be published without reference to any derivation method (as a 
calculational trick that just worked), which gave me a possi- 
bility to "bootstrap" the publications of the theory of AO — 
a non-trivial task given a complexity and utter novelty of the 
theory (it deviated from the then prevailing BPHZ paradigm 
from the ground up; for more on this see [19]). In the final ver- 
sion of the theory, no reference to the R* -operation is neces- 
sary, of course, and the latter remains just an exotic applica- 
tion. In particular, no knowledge of the R* -operation is re- 
quired to understand the example. 

Unfortunately, I had to rely on the method of dimensional 
regularization (regularization-independent understanding of the 
method was achieved after this Letter went to press). The use 
of dimensional regularization somewhat obscured the underly- 
ing elementary analytical mechanism. 

The same example without dimensional regularization 

A regularization-independent treatment of AO to all orders 
was accomplished between 1984 [14] and 1990 (see the sum- 
mary publications [20]). This allows me to explain the example 
of the Letter without the somewhat confusing use of dimen- 
sional regularization which nullifies certain terms that play an 
important role in the mechanism of factorization; the same 
mechanism works to all orders [20]. 

The starting point is now the following 4-dimensional ana- 
log of the integral in (1): 



F = lim 

A^oc 



d 4 p 



1 



(271 ) 4 (p-Q) 2 (m 2 + p 2 ) 



-zln(A/^ uv ) 



0.1 



where: z is a finite constant which can be chosen independent 
of m (this is important; see below); A is an intermediate upper 
cutoff; /i uv parameterizes the arbitrariness of UV subtraction 

(the subscript is introduced to explicitly distinguish it from the 
factorization parameter). 

My purpose here is to simply exhibit the factorization 
mechanism for power correction in the final answer, not ex- 
plain how such a final answer is arrived at via a deterministic 
procedure of AO (for that, see [15]). So I am going to simply 
present the final answer in parallel with the formulas in the 
main text — an expansion of F including one power term in 
which the dependence on the large parameter, Q, is cleanly 
separated ("factorized") from the dependence on m. 

Since the dependence on /z uv is logarithmic, it is irrele- 
vant whether one considers Q as large or m as a small pa- 
rameter. As usual in the theory of AO, it is more convenient to 
work with small parameters, i.e. expand in m — > . 

Instead of the integral (2), we now write: 



F Q = lim 

A->oc 



dV 



1 



1 



(27t) 4 (p-Q) 2 



-zln(A/^i uv ) 



. d 4 p 
(2n) 4 Up-Q) 2 



0(p 2 <fl 2 ) } \ m 2 
(O-0 2 



0.2 



To understand the structure of this expression, compare it with 
the formal expansion of the m -dependent propagator: 



1 



1 



m l + p L 



- + ... 



0.3 



In 0.2, the UV subtraction is inherited by the integral corre- 
sponding to the first term on the r.h.s. of 0.3, p~ 2 , whereas 
the second term, m 2 p^ , is accompanied by another subtrac- 
tion and the cutoff 8(p 2 >fl) , which play a role similar to that 



F.V. Tkachov, Phys. Lett. 125B (1983) 85 



4 



of dimensional regularization and minimal subtraction in the 
main text, jl will be seen to be the factorization parameter. 
Eq. (5) is now to be replaced by the following expression: 

z(m) = J 



(2»)« 



1 



1 



m z +p z 



+>(p 2 >^ 2 ) 
p 



0.4 



which is finite both at small and large p. One easily finds by a 
rescaling p —> mp that 

z(m) = m 2 [c 1 ln(m/^) + c ], 0.5 

where both c, are numerical constants. 

A similar power-and-log dependence on Q 2 (with an over- 
all Q- 1 ) is obtained for the second integral in 0.2 via rescaling 
P^\Q\ P - 

Now the "factorized" expansion is as follows: 
1 



F=F + f 



^—.^rnmp^o^ 



= F + z(m)xf 



d 4 p 



(2n) 4 {p-Qf 



{S(p)} + o(m 2 ). 



0.6 



(In fact, the remainder is not just o(m 2 ) [i.e. vanishes faster 
than m 2 ] but is actually 0(m 4 lnra) .) 

Correctness of the result 0.6 (i.e. the indicated smallness of 
the remainder) is easily checked by combining the r.h.s. and 
l.h.s. and simple algebraic rearrangements of various terms. 

The r.h.s. of 0.6 together with 0.4 and the r.h.s. of 0.2 rep- 
resents a decomposition of the initial integral 0.1 into pieces 
(without explicit integrations) in such a way as to obtain an 
expansion in powers and logarithms of the expansion parame- 
ter (small m or, equivalently, large Q). Prior to discussing the 
interpretation of various terms in 0.6, let us focus on the pat- 
tern of subtractions in the above formulas. 

"Lazy cutoffs" 

(i) Note how the cutoffs like d(p 2 >jj 2 ) enter the various 
terms: d(p 2 >H 2 ) is only introduced precisely into those sub- 
tracted terms that correspond to the singularities of the formal 
expansion, whereas the non-singular term p~ 2 carries no cut- 
off. For instance, the expression for z(m) should be regarded 
(this interpretation remains valid in all orders of perturbation 
theory) as obtained by subtracting from the integrand all those 
and only those asymptotic terms for p — > oo that are responsi- 
ble for UV divergence, whereas the cutoff is introduced only 
into the term corresponding to the logarithmic singularity 
which necessarily also has a logarithmic singularity at p = , 
is why the cutoff is needed in the first place. In other words, 
the finiteness of the corresponding integrals is ensured by sub- 
traction of power-behaved terms rather than a hard cutoff. This 
mechanism can be conveniently called lazy cutoff . 

(ii) Recall how factorization is discussed in the language of 
"factorization theorems" (cf. [21], [22]). There, one discourses 
(cf. [21]) upon contributions of various integration regions, and 
the factorization parameter jl is interpreted as the boundary of 
various integration regions. However, a precise mechanism of 
factorization as exhibited in our example 6 is such that the 
"factorized" integrals collect contributions from all p. 

\ (iii) So, the interpretation of factorization in terms of ; 
; "regions" is fundamentally flawed in one trivial respect: infor- j 



; mation in an integral comes not only from the region but first 
; and foremost from the integrand. This flaw manifests itself at 
i the level of non-leading powers: 



(iv) Indeed, if one starts with splitting the initial integral into 
two regions: \p\<\/i\ and lpl>l^il (recall that in this example 

we work with a Euclidean integral), and performs the corre- 
sponding expansions ( Q — » oo and m — » , respectively) then 
the results would contain various (negative and positive) pow- 
ers of [l, and keeping track of their eventual cancellations as 
well as reassembling the correct factorization result 0.6, is hard 
already in one loop and is a hopeless endeavor in higher orders 
of PT where the overall mechanism becomes rather delicate; 
cf. [20]. 



(v) Plainly put, the mechanism of factorization based on hard 
cutoffs fails with power-suppressed corrections and — unlike 
an implicit belief behind the "theory of renormalons" [21] — 
cannot be a basis for a meaningful theory of such corrections. 



Discussion of Eq.0.6 



6 Expansions in perfectly factorized form are essentially unique [15], so 
one cannot talk about different equivalent representations here. 



(i) Eq.0.6 is in exact correspondence with Eq.(6): Tl(Q,0) 
is exactly the first integral in the expression for F , the r.h.s. 
of 0.2; m 2 C\(Q,ll)Qr 2 corresponds to the second integral on 
the r.h.s. of 0.2; C2(Q,/J.)Q~ 2 is exactly the integral with S- 
function on the r.h.s. of 0.6; and the aggregate (i?r<I> 2 (0)) is 

identified with the z(m) , Eq. 0.4. 

(ii) The integral with the 5-function exactly corresponds to 
the much used formulas for coefficient functions of OPE found 
in a companion Letter [6] and worked out in [9] (cf. the 3-loop 
calculation [23] that defined the state of the art in DIS OPE 
calculations). 

(iii) The correspondence {RT<i> 2 (0)) <-> z(m) is a special 

case of the general results on the connection of diagram-wise 
UV renormalization (the i?-operation) and the subtraction of 
asymptotics at large loop momenta (properly interpreted) di- 
rectly from non-integrated integrands [24], [20]. 

(iv) It is a characteristic feature of the method of AO that the 
integrals one obtains in the "factorization" formulas allows a 
direct diagrammatic interpretation in terms of some operators 
— one never needs any matching conditions. 

(v) The all-orders generalization of the discussed example as 
described in the main text after (6) was achieved in [14]— [1 8] 
where, moreover, the results were extended to arbitrary models 
and arbitrary Euclidean asymptotic regimes, as promised in 
remark (iv) in sec. 3 of the main text. 

(vi) An important update to the remark (i) in sec. 3 of the 
main text is that the dimensional regularization is not crucial 
for achieving the perfect factorization: the analysis of [20] 
shows how the required subtractions (with lazy cutoffs) are 
rather naturally implemented in a regularization-independent 
manner. 

(vii) In cases such as gluon condensates, the masslessness of 
perturbative gluons prevents non-zero condensates from ap- 
pearance in lowest orders (due to absence of a mass scale). 
However, non-zero contributions re-emerge in higher orders 
(massive quark self-energy insertions). 

(viii) A theoretical estimate of the aggregates such as vacuum 
condensates comprising the large logarithms in power- 
suppressed corrections requires non-perturbative methods. For 
Euclidean type objects (e.g. local vacuum condensates), the 
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lattice QCD may provide such estimates. In non-Euclidean 
situations, no meaningful method exists. 

Scenario for non-Euclidean regimes 

The scheme used in the main text of the Letter can be sum- 
marized as follows: 

— perform "perfect factorization" of small parameters such as 
quark masses; 

— identify contributions which give rise to large non- 
perturbative logarithms; 

— treat such PT-uncalculable contributions as phenomenologi- 
cal parameters. 

In other words, the small masses etc. are used as a test of sta- 
bility of perturbation theory in order to identify the structure of 
potentially large non-perturbative corrections. 

The argument that the condensates one thus obtains are not 
"the" condensates of [1] is beyond the point for two reasons: 

1) the large lnm contributions cannot be phenomenologically 
distinguished from contributions of, say, instantons; 

2) the usual analyses such as [1] also rely on PT results in re- 
gard of the structure of power-suppressed corrections (the 
OPE used in [1] was the perturbative OPE somewhat dramati- 
cally misinterpreted due to ignoring the non-zero large lnm 
contributions which in Zimmermann's version of OPE are hid- 
den in coefficient functions); in this respect our scenario is only 
a straightforward clarification of the reasoning behind [1]. 

This scenario is applied to non-Euclidean situations in a 
straightforward manner (as was mentioned in [26] in connec- 
tion with the physics of jets). All one needs is a generalization 
of the method of AO that allows treatment of arbitrary asymp- 
totic regimes and phase space integrals along with loops; such 
a generalization was found in [25]. One proceeds as follows: 

• One starts from complete diagrams with non-zero masses 
for quarks, . . . , non-zero p ± , etc. 

• One performs diagram-by-diagram factorization of the small 
parameters (say, < 1 GeV ) using the algorithms of non- 
Euclidean AO [25]. Note that the method of AO yields 
"perfect" expansions in which the large logarithms of small pa- 
rameters are localized in certain integrals. This is in contrast 
with the conventional methods; e.g. the version of OPE due to 
Zimmermann had such logarithms hidden in coefficient func- 
tions, which phenomenologically is incorrect (see [19]). 

• Among various pieces of the resulting factorized expression, 
one identifies (using scalings like those used to obtain 0.5) the 
integrals which comprise the large logarithms of small masses. 
The simple argument lng/m— >ln<2/ J u-lnm/ i u shows that if 

one normalizes a s at jJL of order Q, a typical scale of the proc- 
ess being studied, then the logarithms of masses are automati- 
cally large, so neither analysis, nor even knowledge of the cor- 
responding evolution equations is really necessary (the data on 
power corrections are never precise enough to warrant that). 

• One classifies independent "non-perturbative" integrals (i.e. 
ones that give rise to large logarithms) which would then cor- 
respond to independent operator objects (being lower-orders 
contributions to matrix elements of the latter). Note that inter- 
preting such objects in, say, coordinate representation as local 
etc. operators serves no useful purpose except, perhaps, as an 
inessential notational convenience. Anyway, even if one de- 
cides that one needs such an operator interpretation, one only 
needs to recognize the emerging integrals (by simply looking at 
them) as corresponding to some or other operator objects {not 
necessarily local operators) — one does not need to massage 



them in any way. Their regrouping in a gauge invariant way 
may be necessary though (e.g. via studying dependence on the 
gauge parameter). 

• Note that for identifying the contributing operators, match- 
ing conditions are never needed. Note also that if one follows 
the prescriptions of [25] carefully, one never has to invent a 
priori prescriptions for renormalization of the resulting non- 
local operators (cf. [32]) — the pattern of necessary subtrac- 
tions is dictated by the method of AO unambiguously (apart 
from inessential finite ambiguities similar to the difference 
between MS and MS-bar renormalization schemes). 

• Upon identification and classifications of independent "non- 
perturbative" integrals, one simply treats them as phenomenol- 
ogical constants. 

• One should replace such integrals with constants as a 
whole, i.e. without attempting to scale out powers of small 
masses because such a possibility is a special feature of pertur- 
bation theory. That a typical scale in such integrals should be A 
rather than, say, quark masses was clear since at least the 
studies of QCD via Dyson equations in the early 80s [29]. Arti- 
ficial arguments such as based on Borel resummation of Borel- 
non-resummable QCD PT series [21] cannot add anything new 
here as a matter of principle. 

• If one has trouble identifying independent integrals, the best 
solution is to consider the next order of PT (one does not need 
to do the integrals explicitly, only write out the prescriptions of 
AO [25] explicitly, which requires only patience) where radia- 
tive corrections would allow one to distinguish different op- 
erator objects. 

• To facilitate investigation of contributions involving gluon 
operators (which to lowest order may be zero due to perturba- 
tive masslessness of the gluon), one may force an early appear- 
ance of non-zero gluon operators by introducing a gluon mass 
(concerning gauge invariance of such a procedure consult [27]). 
This is just a trick because quark self-energies would give non- 
zero contributions from higher orders anyway. 

• The examples I studied show that the breakdown of dimen- 
sional regularization observed in [28] in a semi-artificial ex- 
ample, is quite typical in the case of power corrections to, say, 
jet observables, so one must exercise care with intermediate 
regularizations. 

• The success of the QCD sum rules method [1] where one 
dealt with essentially Euclidean regimes and was able to 
parameterize many observables with a few phenomenological 
constants, is not likely to be repeated in non-Euclidean situa- 
tions where one encounters a wider variety of non-perturbative 
operator objects in addition to local operators and has to deal 
with unknown phenomenological functions of kinematic argu- 
ments, sometimes in infinite number. 

• Moreover, whereas vacuum condensates (that are the only 
type of non-perturbative coefficients in Euclidean situations) 
can be estimated from, say, the lattice QCD, for the non- 
Euclidean situations no useful and meaningful method exists. 

• Finally, a critical examination of renormalons-based results 
demonstrates that their concrete content is always traced back 
to (rather imperfect) analysis of diagrams, whereas the renor- 
malon phraseology effectively reduces to the unoriginal propo- 
sition that the corresponding aggregates are non-perturbative 
without providing any clue as to how one could compute them. 



